The Chalker-Coddington Network Model is Quantum Critical 
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We show that the localization transition in the integer quantum Hall effect as described by the 
Chalker-Coddington network model is quantum critical. We first map the anisotropic network model 
to the problem of diagonalizing a one-dimensional non-Hermitian non-compact supersymmetric lat- 
tice Hamiltonian of interacting bosons and fermions. Its behavior is investigated numerically using 
the density matrix renormalization group method, and critical behavior is found at the plateau 
transition. This result is confirmed by a generalization of the Lieb-Schultz-Mattis theorem. 



PACS numbers: 73.40.Hm, 71.30. +h, 72.16.Rn, 75.10.Jm 



Transitions between plateaus in the integer quantum 
Hall effect provide the clearest example of quantum- 
critical behavior in a disordered system. Understand- 
ing such critical points is a challenging problem because 
fluctuations occur over many decades of length and time 
scales, and averages over different realizations of the dis- 
order must be carried out. Critical behavior was pre- 
dicted by Levine, Libby, and Pruisken Q and was ob- 
served experimentally by Wei et al. S for temperatures 
close to the critical point at absolute zero. Progress to- 
ward a theoretical understanding of the plateau transi- 
tion was achieved with the introduction of a quantum 
tunneling network model by Chalker and Coddington 
||. Subsequent numerical studies Q of the Chalker- 
Coddington model yielded values for the correlation 
length exponent v w 2.3 which were consistent with ex- 
periments. To the best of our knowledge, however, there 
has been no exact proof that the Chalker-Coddington 
model is quantum critical. 

The method of supersymmetry (SUSY) can be used to 
analytically carry out disorder averages We apply 

the method to the anisotropic network model and then 
use the Density Matrix Renormalization Group (DMRG) 
algorithm jl0| to study the resulting spin chain. Unlike 
usual spin chains such as the spin-1/2 Heisenberg antifer- 
romagnet, the on-site Hilbert space in the SUSY chain is 
infinite-dimensional. To access the critical point numer- 
ically therefore requires a double-extrapolation to large 
on-site Hilbert spaces and to large chain lengths. We 
use our numerical results to motivate a generalization of 
the theorem of Lieb, Schultz and Mattis (LSM) which 
confirms quantum criticality. 

1. SUSY Spin-Chain for the Network Model. The 
anisotropic Chalker-Coddington model can be repre- 
sented Jll],[l2| by an independent-particle Hamiltonian 
which describes a chain of L (even) edge states alternat- 
ing in propagation forward and backward in imaginary 
time t. Random complex tunneling amplitudes tj(r) be- 
tween adjacent edge states account for the Aharonov- 
Bohm phases accumulated by the electrons as they cir- 
culate around equipotential contours of the random po- 



tential: 



H ip = fdr J ^(-1)^)( T ) id T <k(r) 
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The (—1)-' term has its origin in the alternating propa- 
gation of adjacent edge states and this factor reappears 
several times in the equations which follow. The disorder 
average of the tunneling amplitudes is given by 



t*(r) t f (r') = Jj ho' 5{t-t') 
J,- = [1 + (-l)'JZ] . 



(2) 



The relevant dimerization parameter R = ±1 deep in- 
side the plateaus; the transition between the plateaus 
occurs at R = 0. Disorder averaging of the correspond- 
ing functional integral is made possible with the use of 
SUSY [^3 14 1 as the partition function Z = 1 for each re- 
alization of the disorder. Transfer matrix formalism can 
be used to resolve normal-ordering ambiguities [jE5 16 
and the resulting effective SUSY Hamiltonian may then 
be extracted. It describes interacting spin-up and spin- 
down fermions c a and bosons b a (two spin species are 
introduced to permit the calculation of the disorder- 
averaged product of retarded and advanced Green's func- 
tions which determines the conductivity): 
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Parameter T) > ensures convergence of the non-compact 
bosonic sector and defines the advanced and retarded 
propagators; also the signs g a are given by: 



9a 



1; o = 1,2,10,12,14,16 
-1; o = 3,. ..,9,11,13, 15 
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In Eq. | we have introduced 16 spin operators, the com- 
ponents of a 4 x 4 superspin matrix: 
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Bosonic- valued operators S* 1 , . . . , S s make up the sym- 
metric sector of the Hamiltonian while fcrmion-valued 
operators S 9 ,...,S le are in the antisymmetric sector. 
Despite the fact that H is non-Hermitian, it only has 
real- valued eigenvalues. 

The Hamiltonian commutes with four (fermion- valued) 
supersymmetry generators, [H, Qi a ] = [H, Q20] = 0, 
where 



3 

q 2ct =e[(-i)^w+4av 



(6) 



It is not difficult to see that the supersymmetric Hamilto- 
nian must have a unique, zero-energy, ground state. The 
right and left (ground) eigenstates are therefore annihi- 
lated by the Hamiltonian: H\^ ) = {^ \H = 0. Also, 
the ground state is annihilated by the SUSY charges: 
Qitrl^o) = ( ?2(t|^ , o} = 0. All excited states appear in 
quartets or larger multiples of 4, half with odd total 
fermion content, and these cancel out in the partition 
function by virtue of the supertrace: 
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where N c is the total number of fermions. In the limit 
of 77 — > 00, |^o) — > |0). For finite r] > 0, however, it is 
a remarkable fact that the normalized ground state is a 
superposition of the vacuum state with unit amplitude 
and an infinite number of zero-norm many-body states 
I J) with differing total number of particles : 

00 

l*o) = |o> + <>Av) \ J); (i\J) = V I, J > 0. (8) 
j=i 



This result can be verified directly by observing that 
|*o) = liiri/3-> oo e~P H |0) and application of powers of H 
to the vacuum state yields only zero-norm states. We 
have also checked numerically, for finite length chains, 
that the vacuum state has unit amplitude when the 
Hilbert space is truncated in a way which respect super- 
symmetry (see below) . It is useful to contrast the compli- 
cated ground state of the non-Hermitian SUSY Hamilto- 
nian with the ground state of the Hermitian SUSY ferro- 
magnet which describes a chiral metal with all edge states 
propagating in the same direction ||. Backscattering is 
absent, localization cannot occur, and the ground state 
of the SUSY ferromagnet is simply the vacuum state. 

2. DMRG Analysis. To simply show that the density 
of states (DOS) is non-vanishing it suffices to remove 
one of the spin sectors, for example the j-spins. The re- 
maining f-spin degrees of freedom are then compact, the 
ground state is the vacuum state |0), the DOS is pro- 
portional to (S 1 ) = —(S 5 ) = 1/2, and the Hamiltonian 
can be exactly diagonalized juj . When both spin sectors 
are included the theory is non-compact and highly non- 
trivial. To make further progress we employ the infinite- 
size DMRG method |l(J. The Hilbert space is first con- 
structed systematically on each site by repeated action of 
the double-creation operator S 3 = Introducing the 

integer level index n = 0, 1, 2, ... we add to the vacuum 
state |0) a tower of states built out of the quartets: 
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Truncations with D = An + 1 states preserve supersym- 
metry, as the SUSY generators, Eq. ^, intermix the quar- 
tet of states, Eq. M, separately within each level of the 
tower without changing the total number of particles. 
The DOS remains unchanged, and the ground state en- 
ergy is exactly zero, providing a valuable check on the 
accuracy of the DMRG algorithm which incurs errors 
when, as the chain length increases, the Hilbert spaces 
of the blocks grow beyond the finite limit of M states. 
Increasing M up to limits set by machine memory and 
speed yields systematic improvement in the accuracy of 
the DMRG algorithm. In results reported below we have 
checked that M is sufficiently large to ensure adequate 
accuracy; even in the challenging case of D = 13 and 
M = 170 the ground state, when targeted, had an en- 
ergy which deviated from zero by less than 0.003. 

Reduced density matrices for the two augmented 
blocks, each of Hilbert space size D x M, are formed 
by computing a partial trace over half the chain. For the 
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left half of the chain the density matrix is given by: 



DM 

Here = (i,i'|\&) are the real-valued matrix elements 
of the targeted many-body wavefunction projected onto a 
basis of states labeled by unprimed Roman index i which 
covers the left half of the chain and primed index i' which 
covers the right half. To compute the gap, ^ is chosen to 
be one of the lowest-lying right eigenstates of H . All of 
the eigenvalues of p are real and positive; these are inter- 
preted as probabilities and the (D — 1)M least probable 
states are thrown away. The desired D —> oo limit can 
be reached by extrapolation, Fig. [[]. Fig. || shows that 
the gap is non-zero in the thermodynamic limit for finite 
fixed D — An + 1 (solid lines) , but as shown in Fig. |l| it 
approaches zero in the D — > oo limit as expected. 
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FIG. 1. Gaps to the lowest excited state for chains with 
open boundary conditions and R = 0. The gap for the 
two-site problem was obtained by diagonalizing Eq. 61 of 
Ref. 16. A straight line is fit to the L — 4 points (with 
M — D). Also shown are gaps at L — > oo which are obtained 
from the extrapolations presented in Fig. 2 for the supersym- 
metric truncations D — 5, 9, 13, and 17. A straight line is fit 
to the last three points. 

Also of interest are non-supersymmetric truncations 
D = An + 2, with the state |4n + 1) selected as the fi- 
nal state at the top of the tower. The special case of 
no bosons, D = 2, with on-site states {|0), c|cj|0)} cor- 
responds to the ordinary spin-1/2 Heisenberg antiferro- 
magnet as can be verified by making separate particle- 
hole transformations on the even and odd sublattices: 
C2j-\ <-» Gyt and C2j+ix <-> cJy^,; consequently the gap 
vanishes in this case. As seen in Fig. || the gap also 
vanishes for D = 6. We prove below that all truncations 
with D = An + 2 are gapless. 




a D=2, M=128 
o D=6, M=192 
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FIG. 2. Gaps, for r\ = 10 -4 , extrapolated to L — > oo. 
Truncations which respect supersymmetry (solid lines) and 
non-supersymmetric truncations (dashed lines) are plotted. 
Points at L = oo are obtained by fitting the gap to the form 
a + b/L + c/L 2 . 

3. Licb-Schultz-Mattis Theorem. For half-odd-integer 
spin antifcrromagnets on a periodic chain of length L 
sites (St ee S ) LSM showed that cither (1) the 
ground state is degenerate or (2) there are gapless spin 
excitations in the thermodynamic limit L — > oo. LSM 
employed a variational argument by introducing the uni- 
tary slow-twist operator, 

l/EEexpj^ jr jS A, (11) 

which has the property that W[H,U] = 0(1/ L). Now 
(^ol^l^o) = because U — » — U under parity of reflec- 
tion about the middle site (j *-* L — j) combined with a 
rotation of 180 degrees about the y-axis (S z <-» — S z ). So 
U either creates a low-energy excitation above a unique 
ground state or mixes degenerate ground states. In con- 
trast, for integer spins U — » U under parity and U docs 
not, in general, create a low-energy excitation or switch 
degenerate ground states. 

In the SUSY problem we are able to make a stronger 
statement, because we know that the ground state is 
unique by supersymmetry. The natural generalization 
of the LSM slow-twist operator for the SUSY chain is: 

U = exp j J pj(-l) j [n c (j) +n b (j) - 1}\ , (12) 

where rib(j) = b^ a bj a is the number of bosons on site j 

and n c (j) ee c^Cja is the number of fermions. It reduces, 
in the D = 2 limit of no bosons, and after the particle- 
hole transformation is taken, to the usual LSM opera- 
tor Eq. |ll|. In the ground state, the sum n c (j) + ni,(j) 
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is always an even number, so U respects the periodic 
boundary condition and it also has the desired property 
U'[H,U] = 0(1/ L). The canonical parity transforma- 
tion now takes the form: 

3 L-j 

(-iy <-> (-iy 

b 3 -\ <-> b L -ji (13) 

again reducing in the absence of bosons to the usual LSM 
parity operation. The supersymmetric Hamiltonian is 
invariant under this operation, but only at the presumed 
critical point R = r\ = 0. The slow-twist operator U, 
however, changes form because while n c — > 2 — n c , the 
number of bosons rib remains invariant: 

0---exp J^{-iy[n c (j)-n b (j)-l]^. (14) 

However, U is invariant under global supersymmetry ro- 
tations, [Qio, U] = [Q2a, U) = 0. As only the ground 
state is a SUSY singlet, U cannot create a low-energy 
excitation for truncations which respect SUSY, consis- 
tent with the above DMRG results. Indeed, from Eq. || 
it follows that | (* |C7|* > I = 1 and thus U\^ ) does not 
contain a component orthogonal to the ground state |\&o)- 
(This can be viewed as an alternative proof of Eq. pi) For 
non-SUSY truncations D = 4n + 2, however, the ground 
state does not obey Eq. || instead |(\E'o|k r | , I , o}| < 1 as can 
be verified either directly for small chains, and in the spe- 
cial case D = 2 of no bosons (for which (^ol^l^o) = 0), 
or by appealing to the fact that U changes form under 
the parity operation, Eq. |lj. For sufficiently large D 
and r\ > the ground state approaches the SUSY ground 
state and is thus unique; therefore U creates low-energy 
excitations Jl9[| , and chains with the non-SUSY trunca- 
tion D = An + 2 are gapless in the thermodynamic limit. 

Before examining the physically relevant D — > oo limit, 
first consider the large spin limit of one-dimensional 
nearest- neighbor Heisenberg antiferromagnets. Chains 
with spin S = n+1/2 (and even- numbered D = 2S + 1 = 
2n + 2) are gapless for all integer n > 0. Chains with 
S = n + 1 (and odd-numbered D — 2n + 3) are Haldane 
gapped, but this gap must vanish in the limit of large spin 
to accord with the gaplessness of the half-odd-integer in 
the n — ► oo limit. This reasoning can be checked by 
a simple rcnormalization-group argument using the beta 
function for the non-linear a-model which shows that the 
gap for integer-spin chains vanishes as e~ . Likewise, 
for rj = + , continuity requires that SUSY truncations 
with odd- numbered D = 4n+l must converge to the gap- 
less behavior exhibited by the non-SUSY even-numbered 
D = 4n + 2 truncations in the n — > oo limit. Hence the 
Chalker-Coddington model is quantum critical at R = 0. 
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